Stochastic Formation of Polariton Condensates in Two Degenerate Orbital States 
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We explore the exciton-polariton condensation in the two degenerate orbital states. In the honey- 
comb lattice potential, at the third band we have two degenerate vortex-antivortex lattice states at 
the inequivalent K and K'-points. We have observed energetically degenerate condensates within the 
linewidth ~ 0.3 meV, and directly measured the vortex-antivortex lattice phase order of the order 
parameter. We have also observed the intensity anticorrelation between polariton condensates at 
the K- and K'-points. We relate this intensity anticorrelation to the dynamical feature of polariton 
condensates induced by the stochastic relaxation from the common particle reservoir. 

PACS numbers: 



Microcavity exciton-polaritons are a solid state sys- 
tem where bosonic condensation [1] has been explored. 
A semiconductor microcavity with embedded quantum 
wells (QWs) realizes the strong coupling between QW 
excitons and microcavity photons, which results in the 
new eigenmodes of the system called upper and lower po- 
lar itons [2|, [3[. At low densities, polar itons are regarded 
as bosonic particles, and lower polaritons (LPs) are ex- 
pected to exhibit dynamic condensation \4\. LP conden- 
sation has been observed in inorganic and organic 
Q semiconductors. One unique feature of polariton con- 
densates comes from their non-equilibrium nature. Be- 
cause polaritons leak from the cavity as emitted photons 
before reaching the thermal equilibrium, the stochastic 
relaxation plays a role for the formation of the conden- 
sate. Recently the stochastic polarization build-up of the 
polariton condensation has been theoretically and exper- 
imentally explored Io|, HI- 



Using the thin metal deposition method, single trap 
12|, one-dimensional array [13| and 2D square lattice 
potentials [14| have been studied, where the dynamical 
and bottleneck condensation in meta-stable higher bands 
has been observed. A band gap between the bands sup- 
presses a LP relaxation to the lower band, which results 
in the meta-stable condensation at the top of the gap. 
Due to the short lifetime of the LPs, they can escape 
from the cavity as emitted photons prior to the relax- 
ation to the lower bands, which allows the direct ac- 
cess to the particle density and phase distribution of the 
higher band condensates. The mode competition and 
the relaxation dynamics for different energy modes have 
been studied experimentally and simulated by the rate 
equations [HI, (3- For the degenerate energy states, the 



relaxation to these states are initiated by the indepen- 
dent spontaneous scattering and the stochastic popula- 
tion build-up will be induced between these states. Here 
we report the stochastic formation of polariton conden- 
sates in two degenerate orbital states in the honeycomb 
lattice potentials. 

The honeycomb lattice potential is implemented by the 
thin metal deposition method with ~ 150 jaeV potential 
strength and the period a ~ 2 /im (Fig. QJa)). The 



characteristic kinetic energy is 



h z k 
2m v 



1.5 meV with the 



polariton mass m p: the reduced Plank constant H and 
the unit wavenumber fen = -77? — at a red detuned area 
(A - -2 meV). The Brillouin Zones (BZs) with three 
high symmetry points (r, M and K) are shown in Fig. 
DJb). The band structure (Fig. [He)) is calculated with 
single-particle plane wave bases for the weak potential 
regime. There are two inequivalent K- and K'- points at 
the vertices of the first BZ with 3-fold rotational symme- 
try. Under the weak periodic potential, this 3-fold degen- 
eracy is split into lower degenerate doublets and upper 
non-degenerate singlet with an energy gap (~ 35 fieV 
with our potential), which suppresses the LP relaxation 
from the upper singlet to the lower doublet states. 

The details of our sample with 12 GaAs QWs in Al- 
GaAs/AlAs planer microcavity are shown in the previous 
paper Q . All measurements in this paper are performed 
at ~ 4 K, and around a red-detuned area with the detun- 
ing A ~ -2 meV. LPs are created by Ti:Sapphier laser 
of ~ 767.5 nm (QW exciton resonance) incident at ~ 60 
degree in the pulsed mode with a 76 MHz repetition rate 
and ~ 3 ps pulse width. 

We first explore the far-field polariton distribution. 
Sharp intensity peaks in the far-field pattern shown in 
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FIG. 1: (a) A photograph of the sample surface. Brighter 
areas where the thin metal films are deposited correspond to 
the high potential regions, while darker areas with no metal 
correspond to the polariton traps with low potential, (b) Cor- 
responding Brillouin Zones (up to the 3rd) and high symme- 
try points, (c) Calculated band structures with the potential 
amplitude V = 0.15£ , where E = ti 2 (2tt) 2 /2ma 2 . 



Fig. [2fa) indicate that large LP population condenses 
above the threshold pump power P/Pth ~ 2 at the K- 
and K'-points in the 3rd BZ. We note that LPs are more 
or less equally populated at all six points within ~15% 
difference. At the K-points, the lowest three ort honor mal 
eigenstates are split into two degenerate doublet ) 
and 1^2") an d the upper singlet l^f") [l5[. The pump 
power dependence of the K-point energy, linewidth and 
intensity are shown in Fig. [2fb). Across the threshold 
pump power (P t h ~ 4 mW), the intensity increases non- 
linearly. The energy is blue-shifted at the threshold, and 
gradually increases with increasing the pump power. The 
linewidth shrinks at the threshold and increases above 
the threshold due to the polariton-polariton interaction. 
All of above results show the characteristic behavior of 
polariton condensation 0, 0, [l6| . It is also experimentally 
confirmed that the energy of all six K- and K'-points are 
degenerate within their spectral linewidth (~ 0.3 meV). 
The energy-resolved far-field patterns reveal the relax- 
ation path of LPs (Figs. E(c)-(e)). From the 2nd T- 
point to the K(K / )-points, LPs relax through the 3rd BZ 
not the 2nd BZ, which is consistent with the condensa- 
tion at the 3rd band K(K / )-points. Because LPs relax 
by polariton-polariton interaction and polariton-acoustic 
phonon interaction in low temperatures both of which 
conserve the in-plane momentum, it is natural to think 
that LPs relax to either or ) state (not super- 



position states of two). 

The order parameter of l^f") and ) condensates 
are directly measured by a modified Mach-Zehnder inter- 



FIG. 2: (a) A time-integrated far- field image (in linear scale) 
showing LP condensation in momentum space above the 
threshold pump power; P/Pth ~ 2. The white lines indicate 
the BZs, which are calculated with the diffraction signals of 
the reflected laser, (b) The pump power dependence of the LP 
population, the energy shift at k=0 measured from the value 
below the threshold P = 2.5 mW and the linewidth of the K- 
point condensate. All the values are averaged over the six K 
and K'-points, and the errorbars show the minimum and the 
maximum value of them, (c) (d) (e) Energy-resolved time- 
integrated far- field images at 1612.06 meV (c), at 1611.88 
meV (d) and at the K-point condensate energy 1610.49 meV 
(e) taken at P/P th ~ 2. 



ferometer [18|. In a signal arm, signals from all six K- 
and K'-points are chosen. In a reference arm, we choose 
one of the three K- or K'-points using a 150 jam pin- 
hole (corresponding to Ak ~ 0.31 /im" 1 ) as a reference 
plane wave for the interferogram. We use a band-pass 
filter with 1 nm pass band width to filter out the resid- 
ual light from other states as well as the scattered laser. 
The near-field image observed in the signal arm shown 
in Fig. 09(a) is consistent with the theoretical particle 
density distribution of ) and ') in Fig. [3fb) with 
the intensity peak-to-peak distance of 2>/3 /im. We note 
that 1^3") and ) have the same intensity distribu- 



tion. This result shows that 



and ) are inco- 



herent in the time-averaged sense, because superposition 
states have the different intensity distribution [19]. The 
theoretical phase profiles of l^f") and \^f^ ) are shown 
in Fig. 09(c) and Fig. [3f d) respectively, where clockwise 
and counter-clockwise phase rotation exist to form the 
honeycomb lattice geometry. At the position where the 
phase rotation exists the intensity goes to zero (in Fig. 
09(b)). This vortex- ant ivortex lattice order is directly de- 
tected in the experimental interferogram shown in Figs. 
09(e) and (f). By extracting the off-axis components of the 
interferogram [2o|, HH , we can reconstruct the phase pro- 
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FIG. 3: The time-integrated interferogram for the K- and K'-points condensate at P/Pth ~ 2. (a) The near- field image which 
includes all the K and K'-points signal, (b) The theoretical near-field image of ) and ). (c) (d) Theoretical phase 
profile of ) (c) and l^f" ) (d). (e) (g) The experimental interferogram between the signal from all the K and K'-points in 
the signal arm, and the signal from only one K-point in the reference arm (e) and the extracted phase profile of its off-axis 
component (g). (f) (h) The same interferograms, but for the K'-points. 



files of the \^^) an d ) order parameter (Figs. [3f g) 
and (h)), which well match with the theoretical phase 
profile shown in Figs. [3fc) and (d). By selecting the 
K-point or K'-point reference, we can extract the order 
parameter profile of l^ff) or ) condensates respec- 
tively. This vortex- ant ivortex lattice order is only seen 
above the threshold, and thus originating from the K- 
and K'-points condensates (22| . 

Finally we explore the dynamics of the mode compe- 
tition between the degenerate l^f") and l^f" ) conden- 
sates. Figure SJa) shows the normalized second order 
auto and cross correlation function of l^f") and l^f- ) 
condensates with the pump power dependence of the T- 
point and K(and K')-points condensate intensity (Fig. 
Egc)). #( 2 )(K,K,t = 0) (^ 2 )(K',K',r = 0)) is the mea- 
sured auto-correlation function of one K(K / )-point se- 
lected by the pinhole, while g( 2 )(K,K',r = 0) is the 
measured cross correlation function between one K-point 
and one K'-point. In both cases, the band pass filter 
with 1 nm band width is used to filter out the resid- 
ual signals and Si avalanche photodiodes are used for 
photon-counting. We can see the anti-bunching behav- 
ior between l^^) and l^f- ) condensates. The degree 
of the anti-bunching has the maximum at P/Pth ~ 2, 
where the K(and K'j-points condensates are dominant 
in the system. With increasing the pump power, the 
T-point signal becomes dominant and the anti-bunching 
becomes smaller. We relate this intensity anticorrelation 
to the mode competition resulting from the stochastic 
formation of l^f") and ) condensates. In each pulse, 
the l^f") and ) condensates have the common par- 
ticle reservoir. The LP relaxation is first initiated by 
the spontaneous scattering through the LP-LP and LP- 
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FIG. 4: (a) Experimental results for the normalized second 
order auto correlation, # (2) (K,K,r = 0), # (2) (K',K',t = 0) 
and the cross correlation </ 2 )(K,K',t = 0). (b) Simula- 
tion results for <? (2) (K,K,r = 0), <? (2) (K', K', r = 0) and 
# (2) (K,K',t = 0) as well as # (2) (A11, All, r = 0). (c) (d) Ex- 
perimental (c) and simulation (d) results for the pump power 
dependence of the intensity at T-point and the K-point con- 
densate, (c) The peak intensity values at the T-point and one 
K-point are plotted, (d) The T-point intensity (|^g| 2 ) and 
the K-point total intensity (|^k| 2 + |^k'| 2 ) are plotted. 



phonon scattering. When the population at one state 
first reaches the quantum threshold, the stimulated scat- 
tering to that state is turned on and the relaxation into 
that state is more enhanced than the other, while the 
total particle number is determined by the pump power 
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above the threshold. Thus, the intensity anticorrelation 
between the l^f") and ) condensates occurs. 

We also simulate this stochastic condensates formation 
using the complex- number Langevin equations, similar to 
those for the polariton condensate order parameter with 
spins [Io|, Hlj . We extend the model to include the T- 
point condensate. We assume that the relaxation from 
the l^f") and ) condensates to the ground state is 
mediated by the thermal reservoir because multi-phonon 
scattering processes are included. The set of equations 
are 



d 
d 

dt* G 



R 



I 2 )] 

(2) 



Jt NR 



\ [-r <T - 7c (|M/ (T | 2 + i) +7 iv. _ 

i [-r G + -f G N R + -y c (\y K \ 2 + \t> K , 
0GR(t) + e GK (t) + e GK ,{t), 
p-r R N R -~fN R (\y K \ 2 + i) 

- jN R (\* K ,\ 2 + l) -TciVfl (|* G | 2 + 1), (3) 

where a = K or K' , ^ k' and are the or- 

der parameter for l^^), ) and the ground state 
(the T-points) condensates respectively [23] . Because our 
pump laser is a pulsed one, we put Nr(0) = J dtP(t) 
and set P=0 for all time [10j. The dynamical thresh- 
old condition Nr(0) = J P t h = ^k/j for the case with- 



out the T-point condensate [24[ is used here, though the 
presence of the T-point condensate may change the ac- 
tual threshold. We numerically solve the equations us- 
ing a (fifth-order) Adams-Bashforth-Moulton predictor- 
corrector method and take 1000 samples for each power 
[iBj . The simulation results are shown in Figs. H^b) and 

(d). 

In the simulation, at P/P t h < 1-9, g^(K,K) and 
q^{K' ', K') show the thermal state behavior, g^ 2 \r = 
0) - 2 [26]. ^( 2 ) (All, All) - 1.1 at P/P th ~ 1.9 indi- 
cates that the order parameters are well defined at the 
K- and K'-points and the states are coherent states. The 
stochastic relaxation induces the anticorrelation between 
tik and uk' ■> and it appears as the bunching in (K, K) 
and g( 2 \K' \K f ) and the anti-bunching in g^ 2 \K,K'). 
g( 2 ) (All, All) ~ 1.1 also shows that there is an excess 
noise induced by the relaxation process from the reser- 
voir. This excess noise makes the degree of the anti- 
bunching in g( 2 \K,K') smaller than the degree of the 
bunching in g^ 2 \K,K) and gW(K',K'). We note that 
the thermal and quantum depletion [lj, which are not 
included in this model, are other important sources for 
the excess noise in experiments [I?} . With increasing the 
pump power, the T-point condensate become stronger 
than the K-point and K'-point condensates (Fig. Iljd)), 
and the excess noise is also increasing (Fig. H^b)). Here 
the relaxation from the l^f") and ) to the T-point is 
comparable to the stimulated scattering from the reser- 



and ) condensates. The small degree of the anti- 
bunching in experiments indicates that not all reservoir 
particles are common for l^f") and l^f- ) . Because of 
the momentum conservation in the relaxation process, 
some LPs depending on their momentum can only be 
scattered to either l^f") or \^/^ ). This simulation qual- 
itatively reproduces the experimental result. 

We briefly note that although the LP-LP interaction 
cannot change the LP population in and ) due 
to the wave-function parity symmetry, the phonon scat- 
tering or impurity scattering can in principle transfer the 
LPs from l^f") to ) and vice versa, and they may 
contribute the intensity anticorrelation. We also note 
that the LP-LP repulsive interaction may induce the co- 
herence between and ) condensates. Such dy- 
namics between the and ) condensates remains 
for the future study. 

Our results show the peculiar characteristics of polari- 
ton condensates in two degenerate orbital states, where 
the stochastic relaxation process plays an important role. 
We also show that it is possible to form the novel phase 
order, which is vortex- antivortex lattice order, in the po- 
lariton condensates originating from the single particle 
wavefunction in honeycomb lattice potentials. 
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voir to |^3 ) and 



i*f>. 



Thus the noise effect becomes 
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bigger and brings the intensity fluctuation to the 



K ) 
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